A continuum model for low temperature relaxation of crystal steps 
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High and low temperature relaxation of crystal steps are 
described in a unified picture, using a continuum model based 
on a modified expression of the step free energy. Results are 
in agreement with experiments and Monte Carlo simulations 
of step fluctuations and monolayer cluster diffusion and re- 
laxation. In an extended model where mass exchange with 
neighboring terraces is allowed, step transparency and a low 
temperature regime for unstable step meandering are found. 



PACS numbers: 05.40.-a, 61 .46.+w, 68.35.Fx, 36.40.Sx 

The upsurge of nanotechnologies and of microscopic 
visualization techniques (such as Scanning Tunelling Mi- 
croscopy) in the past 15 years, raises the challenge for 
modeling crystal surfaces morphology at smaller and 
smaller scales. Surface dynamics crucially depends on 
the underlying microscopic discreteness, especially in 
non-equilibrium or low temperature conditions. Ac- 
counting for the presence of crystal steps has allowed 
some important breakthroughs in the description of 
growth and relaxation of nanostructures [0,^]. In this 
Letter, it is pointed out that steps themselves exhibit 
a discrete sub-structure that drastically affects their 
dynamics at low temperatures, and that a continuum 
model based on a modified free energy allows one to ac- 
count for this low temperature regime. 

An isolated step is always rough: steps do not have 
macroscopic facets or sharp angles (see Ref. Q for a dis- 
cussion). This statement has motivated a modeling for 
mass transport along steps (edge diffusion), first devel- 
oped by Mullins based on a description at length- 
scales larger than the distance between kinks. Measure- 
ment of the time correlations of fluctuating steps Q has 
revealed a low temperature regime for edge diffusion not 
explained by this model. Systematic low temperature 
deviations from its predictions, have also been reported 
in experimental and kinetic Monte Carlo (MC) studies 
of monolayer island diffusion j^^j, and relaxation from 
a deformed shape [0. 

In the following, a continuum model is presented, 
based on a modified expression of the step free energy 
that explicitly accounts for edge-atoms (i.e mobile atoms 
at the steps) . This model exhibits high and low temper- 
ature regimes in agreement with kinetic MC simulations 
and experiments. When mass exchange with neighbor- 
ing terraces is allowed, it also accounts for the recently 
observed step transparency (also called permeability) on 



high temperature Si(lll) surfaces and for low tem- 
perature wavelength selection of unstable step meander- 
ing during growth |d| . 

The free energy of a step is traditionally taken to be 
proportional to its length M : 



dx 7 (l + aC) 2 )) 1/2 



(1) 



where £(x, t) is the meander of the step with respect to 
the straight step configuration, and d x denotes the par- 
tial derivation with respect to x. At low temperatures, 
one can neglect atom detachment from steps to terraces 
|Q. Mass transport then only occurs via mobile edge- 
atom diffusion along the step edge, and is driven by gra- 
dients of the chemical potential \i — ^(S^o/SQ = 
where Q is the atomic area, S/S(^ denotes a functional 
derivative with respect to and 7 = 7 + 7" is the step 
stiffness. Step motion results from the divergence of the 
local mass flux j = —aDT,d x [ijkBl ', where Dr, is the 
macroscopic diffusion constant for mass transport along 
the step. Thus, 



d t ( = d x 



d x (Q 7 k) 



(2) 



Eq.(g) is the usual starting point for studies of edge 
diffusion-driven step dynamics. It is always valid for 
lengthscales larger than the distance between kinks, and 
long enough timescales. 

As an example, let us consider a [110] step on a 
Cu(100) surface, where the kink energy is Ek = 0.13 
eV Q. At low temperatures, kink density is Nk ~ 
2exp(—Ek/kBT)/a, where a is the lattice spacing. The 
distance between kinks is then w 75a for T = 



300K, and NZ 



10 6 a at T = 200K. Hence, not only 



N k S> a, but at low enough T, 



may be much 



larger than observation lengthscales: Scanning Tunelling 
Microscopy nowadays allows one to study step fluctua- 
tions up to atomic scales QJhJ, and monolayer clusters 
of several nanometers are observed Q. Moreover, step 
relaxation timescales related to edge-atom motion from 
kink to kink can also become large when Nk is small. 

The very different role played by mobile edge-atoms 
and atoms incorporated into the solid suggests that the 
step should rather be described as a heterogenous phase 
at scales smaller than NjT 1 . Hence, the total free energy 
shall be written as: 



dx 
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for a step wandering about the closed-packed orienta- 
tion x. c(x, t) is the macroscopic concentration of mo- 
bile edge-atoms along the step, as presented in Fig.l. 
The term (c — c° eq ) 2 accounts for the departure from lo- 
cal equilibrium, and will be seen to be irrelevant at high 
temperature. We shall see in the following how step re- 
laxation based on Eq.(||) accounts for high and low T 
regimes in a wide variety of physical situations. Kinks 
are not explicitly described in this model. Nevertheless, 
the relaxation of c to its equilibrium value c\ q implicitly 
involves the kink distribution, as will be seen later. 
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FIG. 1. The step position £(x,t) and the macroscopic 
edge-atom concentration c(x, t) are both needed in a con- 
tinuum description for low temperature dynamics of crystal 
steps. 

For small perturbations, Eq.(|J) is expanded to second 
order in £ and u = c — c® q . At thermal equilibrium, 
equipartition implies that each Fourier mode Uk or of 
wavevector k carries the same amount of energy: 
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As a first mean field approximation, edge-atoms can be 
considered as non-interacting; thus 



\Uk\ 



^eq t 



and a = ksT/c^. 

Let us now turn to the dynamics. In a way similar to 
model C in critical phenomena |T^| , two evolution equa- 
tions are written. The first one for the non-conserved 
step position £ reads: 



dtC 



= -A 



5JF_ 



v 



(5) 



where A is a kinetic coefficient, and 77 is a Langevin force. 
A second evolution equation is written for the conserved 
total concentration of atoms (i.e. atoms in the solid and 
edge- atoms) -0 = C/^ + c: 



d t ip = d x 



where B is a mobility, and q is a conserved noise. Using 
Eq. (||) in Eqs. a se * °f coupled evolution equations 

is found: 



dtC , 



-eq) 



V 



d t c = d x [Bd x c - q] - u(c - c eq ) - rj 



(7) 
(8) 



where v = Act , and c e 
Thomson relation: 



is found to obey a Gibbs- 



c eq — c eq(l T Tk) . 



(9) 



where k is the step curvature, and T = £l^/k B T is the 
capillary length. The obtained equations share similar- 
ities with the model of Appendix A of Ref. jl3| for line 
diffusion or that of Ref. |l4j for surface diffusion. Nev- 
ertheless, they were not derived from a free energy as 
done in this Letter, and to our knowledge, the results 
mentioned in the following have not been pointed out by 
these authors. Following Ref. the correlations of the 
Langevin forces are found within a local thermodynamic 
equilibrium approximation: 

(77(2:, t)r)(x', t')) = 2vc(x, t) S(x - x')8(t - t') , 

(q(x, t)q(x', t')) = 2Bc(x, t) S(x - x')S(t - t') . (10) 

As a first approach, we now propose some phenomeno- 
logical expressions for the kinetic coefficients. The 
macroscopic diffusion constant of edge-atoms along the 
steps is approximated by that of a tracer edge-atom 
on a frozen step with a given kink density N^. Defin- 
ing the diffusion constant of mobile edge-atoms be- 
tween kinks D and the kinetic attachment lengths d± = 
a(exp(E±/ kgT) — 1), where E± are the additional en- 
ergy barrier (with respect to diffusion) for atoms to stick 
to a kink from both sides, and using the result of Ref. 
ma it is found that: 



B = 



D 



l + N k (d++d-) 



(11) 



The macroscopic attachment coefficient v is the inverse 
of the relaxation time of the concentration, which is re- 
lated to the timescale of diffusion of a mobile atom from 
one kink to another. Hence, 



BNi. 



(12) 



We shall first address the case of a straight step at 
equilibrium fluctuating about the close-packed direction 
x. The time correlation function 



G{t) = ([((x,t)-((x,t + T)] 2 ) 



(13) 



has been measured in experiments and MC simulations 
PJ . This quantity is evaluated within the quasistatic ap- 
proximation, which stipulates that the edge-atom con- 
centration reaches a steady state on time scales much 
shorter than kink motion. Thus the l.h.s. of Eq.(||) van- 
ishes. Comparing the quasistatic and full dispersion re- 
lations as in Ref. , and using Eq. (|l^) , the quasistatic 
approximation is found to be valid when 



Fa 2 - 



C eq « N k 



(14) 
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At low T, in a simple bond-counting model and for 

a step along a high symmetry orientation, Ta 2 c eq ~ 

a 2 exp(2Ek/kBT), where 
fl 1 i„ at. „-i 



a 2 exp(-£; fc /fc s r) « A^ 2 



S fc is the kink energy. At high T, c° 9 ~ 1/a, Afc ~ 
and r — > 0. Hence, the quasistatic limit is valid at all T. 

Linearizing Eqs.(0J|), in the quasistatic limit, G(t) is 
easily evaluated. For large observation timescales, long 
wavelength fluctuations A ^> N^ 1 dominate, and 



Gi on g (r ) 



a 2 r(3/4) 



(15) 



where b 2 = aksT/j is the step diffusivity. This expres- 
sion corresponds to the one given in Ref. J2(]] starting 
from the Mullins model Eq.(||), with Dl = aBc° eq as ex- 
pected from Ref. ||. For short observation timescales, 
only short wavelengths (A <C N^ 1 ) contribute to G, and: 



7 3 / 2 



G short (r) = " („ C °62)VV/ 2 . 



(16) 



Using Eq.(|l2|) and the relation b ~ A& -valid at low T, 
the crossover between the two regimes is found to corre- 
spond to G(t) r~j a 2 . This result was found by Giesen et 
al Q , by means of a discrete random kink model and MC 
simulations. From the relation G(t*) ~ a 2 , the cross-over 
time between the two regimes is found to be: 



t* 



(Nla 2 Bc°y 



(17) 



Using numerical values for Cu(lln) vicinal surfaces given 
in Ref. Q, one finds t* ~ Kr 19 exp(14870/T) s, where 
T is in Kelvins. With observation times t* ~ Is [Q, 
the crossover is found for T « 340K, in quantitative 
agreement with experiments Q. Formally, the cross- 
over from high to low temperatures (Eqs.(|ll^|l6|)) looks 
similar to that from line diffusion to terrace diffusion 
reported by several authors Neverteless, it is 

physically different since there is only line diffusion here. 

Monolayer cluster relaxation and diffusion might also 
be addressed by this model. Let us consider an circular 
island having small perturbations about its mean radius 
Ro defined by R(8) = Rq + p(6). The polar coordinates 
R and 8 are used. In a 'circular model', x is simply 
replaced in the model by Rq9 for small perturbations. In 
the linear approximation, the relaxation time of a small 
perturbation p — e cos(n#) is: 



t n = 2tt 



Rl(u/B) 



vTa 2 c eq n 2 (n 2 



(18) 



A crossover from t r ~ R$ when Ro ^> N^ 1 to t, 
when i?o <C N^ 1 is found in agreement with kinetic MC 
simulations in Ref. 0. 

The diffusion constant of a cluster is defined as: D c = 
( r Cj\/(^))/^' wriere yqm indicates the position of the 



center of mass of the cluster. The cluster diffusion con- 
stant is calculated without using the quasistatic limit, 
and reads: 



Dr. = 



i 



itRq Rl/B + l/v 



(19) 



The high T behavior D c — ^Dl/itRq, calculated from 



Eq.(g) in Ref. 19|l7[], is recovered when R > A fe 



provided once again that D^ = aBc eq . In the low T 
regime, where Rq -C A^ -1 , another scaling limit is found: 
D c ~ i?Q_ x , in agreement with previous experimental M 



or MC H studies. The circular approximation catches 
the essential physical point, which is the existence of 
orientations for which A^ 1 is much larger than the size 
of the cluster. Nevertheless, including anisotropy in our 
model is needed for a quantitative comparison with low 
T experiments and Kinetic MC simulations. 

At higher T or during growth, steps exchange mass 
with terraces. For low kink concentration Nk <C a -1 , 
direct exchange from kink to terrace can be neglected, 
and step meandering is weak, i.e. d x (^ <C 1. An extended 
model may then be written, with Eqs. ([?]), (^) and: 

d t c = d x [Bd x c] ~ v{c -c eq ) + J+ + J_ (20) 
J± = P±C± - v±c (21) 
d t C = D s S7 2 C + F-C/t (22) 

where C is the concentration of adatoms on terraces. 
D s is their diffusion constant, F the incoming flux on 
terraces, and r the adatom desorption time. Langevin 
forces were omitted in Eqs.(p(i|-p2"l) for the sake of clarity 
in this brief exposition. + and — designate the lower and 
the upper sides of the step respectively, v± and (3± are 
kinetic coefficients. At equilibrium c = c\ q and C = C° g , 
and there should be no mass flux (detailed balance) so 
that J+ = J_ = 0. Thus, C° q /c° q = v+/(3+ = v-/(3- 
Since we address the case of weak meandering, exchange 
mass fluxes between steps and terraces are given by 
J± w —D s d z C±, where z is define in Fig. I. These latter 
equations allow to close the model and to evaluate the 
concentration. 

In the case of slow diffusion along steps {i.e. B small), 
or for long wavelength perturbations (larger than A^T 1 ), 
the first term in the r.h.s. of Eq.(p0|) vanishes. If an 
additional approximation is made in taking the varia- 
tions of the adatom concentration C on terraces to be 
much slower than the relaxation time (y + v + + v-)^ 1 
of the edge-atom concentration c, we can set dtc — in 
Eq.(Efjf). The resulting model may be written: 



1 




d t ( = D s d z C+ - D s d z C. 



(23) 



D s d z C ± = ±0±(C ± - C eq ) + f3 a (C+ - C_) (24) 
C eq = C° eq (l + rK) (25) 

where effective kinetic coefficients are defined via 

P+/13+ = 0-//3- = and/? - P+v-/(v+ 
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v+ + v-). With Eqj22|) and Eqs.(|23f|25j), we have ob- 
tained the standard B model for "transparent" steps as 
a special limit. Step transparency (i.e /3q 7^ 0) is un- 
derstood as the possibility for an atom to attach to a 
step from a terrace, and detach to the other terrace be- 
fore reaching a kink. Transparency appears as a natural 
ingredient of the model equations (|?|, pO| - f22| ). 

As a last remark, we shall calculate the most unstable 
wavelength for the meandering instability first addressed 
by Bales and Zangwill |22], with step relaxation provided 
by the full model Eqs.( |20[|2"2] ). We use the following pa- 
rameters: /3_ = and V- = (i.e. no mass exchange 
with the upper terrace); C® q = 0, v+ = 0, and /3 + — > 00 
(this implies C+ = 0). Moreover, desorption is taken to 
be vanishingly small: 1/r = 0. The growth rate of a 
small in-phase perturbation £(x, t) — exp(icj£ + ikx)(^ UJ k 
of all steps on a vicinal surface reads: 

-fc 4 I\Bc° + F(k£t&nh(M) + sech(fcf) - 1) 

= ; — . . , „ (26) 

1 + [B/v)k 2 v ' 

where £ is the mean inter-step distance. The most un- 
stable wavelength is calculated in the long wavelength 
limit kl -C 1 (valid for small fluxes F). One finds: 



Ai = Air(YBc° /Ft 2 ) 1 / 2 



(27) 



when N k 1 <C Ai. In the opposite case N k 1 ^> Ai, a low 
T regime is found, where 



A 2 = 2 1 / 4 7r 1 /2 A }/ 2 7V; 1 



/2 



(28) 



Using activation energies given in Ref . 0] , one finds that 
Ai and A2 follow Arrhenius laws with activation energies 
0.38eV and 0.12eV respectively (within ~ 10% error). 
The low T regime seems to provide the best fit to the 
experimental result of 0.09eV ||. It is not clear though 
how non-equilibrium line diffusion effects pointed out in 
Rcf. jlll combine or compete with these results. 

In conclusion, a model has been presented, based on 
Eq.(||) and Eq.jT^), that accounts both for high and low 
temperature step relaxation dynamics observed in ex- 
periments (4|(| and Kinetic MC simulations . When 
mass exchange with neighboring terraces is added to this 
model, step transparency appears as a natural conse- 
quence of a low kink density. We also point out a low 
temperature regime for step meandering during growth. 

A systematic analysis from a microscopic theory is 
still needed for a more rigorous evaluation of the ki- 
netic coefficients B and v. Moreover, numerical solu- 
tion of the fully anisotropic model is needed in order to 
describe quantitatively monolayer cluster diffusion and 
relaxation. 

Low temperature relaxation of three-dimensional clus- 
ters and nanostructures is a source of long standing con- 
troversies [|. The basic difficulty comes from the singu- 
larity of the free energy for orientations in the vicinity 



of a facet. As opposed to this situation, the free energy 
of a step does not exhibit singularities. Thus, a direct 
generalization of the present study is not possible. Nev- 
ertheless, it provides some milestones for a continuum 
description of three-dimensional clusters |^3| and nanos- 
tructures relaxation. 
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